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Abstract. We review recent literature on the connection between quantum entanglement and cosmology, with an 
emphasis on the context of expanding universes. We discuss recent theoretical results reporting on the production of 
entanglement in quantum fields due to the expansion of the underlying spacetime. We explore how these results are 
affected by the statistics of the field (bosonic or fcrmionic), the type of expansion (dc Sitter or asymptotically stationary), 
and the coupling to spacetime curvature (conformal or minimal). We then consider the extraction of entanglement from a 
quantum field by coupling to local detectors and how this procedure can be used to distinguish curvature from heating by 
their entanglement signature. We review the role played by quantum fluctuations in the early universe in nucleating the 
formation of galaxies and other cosmic structures through their conversion into classical density anisotropics during and 
after inflation. We report on current literature attempting to account for this transition in a rigorous way and discuss the 
importance of entanglement and decoherence in this process. We conclude with some prospects for further theoretical and 
experimental research in this area. These include extensions of current theoretical efforts, possible future observational 
pursuits, and experimental analogues that emulate these cosmic effects in a laboratory setting. 
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The study of the early universe and the large-scale structure of the spacetime has flourished since the introduction of 
the theory of inflation [1] , tested by the observation of anisotropics in the cosmic microwave background [2j [3] . Many 
questions still remain, however. In the late 90s it was found that our universe undergoes an accelerated expansion [4j[5], 
implying that most of the content of the universe is some form of dark energy whose nature remains unknown. Also, 
measurements of the deuterium-hydrogen ratio in interstellar absorption, together with some theoretical models of 
cosmological nucleosynthesis, have provided evidence that ordinary baryonic matter is responsible for only about 20% of 
the matter content of the universe and that there should be some sort of non-baryonic cold dark matter [6] accounting 
for the rest. The abundance of this dark matter cannot be explained exclusively with standard- model neutrinos [7]. 
This mysterious dark matter remains elusive even to the latest progress in experimental particle physics in the Large 
Hadron Collider H Q1 QT] . 

Another important ingredient of a consistent cosmological model is still missing: we do not fully understand the 
nature of the quantum field responsible for the hypothetical inflation period in the very early universe [3]. In the 
absence of a full quantum theory for gravity, which would include backreaction of the quantum fields on the curvature of 
spacetime, quantum field theory in curved spacetime (without backreaction) is the most complete theory so far [12] . The 
gravitational curvature of spacetime has nontrivial effects on quantum fields living on the spacetime when compared with 
their flat-spacctimc counterparts. This is especially interesting in the case of dynamical spacetime backgrounds and/or 
when the spacetime includes horizons preventing an observer from having experimental access to the full quantum state 
of a field |13j . It is known that the gravitational interaction may induce or reveal inaccessible quantum correlations 
in the field state in scenarios such as expanding universes or stelar collapse [21 HH [16] and that particle detectors 
interacting with a quantum field can become entangled |17j in a way that is sensitive to the structure of the spacetime 
that those fields inhabit [T8l[T9] . 

As we will discuss in this paper, the early universe's energy content was dominated by a highly entangled quantum 
field background. If the quantum correlations managed to survive to our days in weakly interacting fields in some form, 
they can provide precise information about the nature and history of the spacetime |20| I21[ 122] . Their study may prove 
useful in constructing models of the early universe and in gaining information about the phenomenon of inflation. In 
particular, we would like to understand the decoherence mechanisms that governed the transition from a highly non- 
classical quantum-correlated field state to the current classical background. By means of this understanding we may be 
able to provide answers to the still unsolved questions about the inflationary theory. 

In this article, we review the different scenarios studied in the literature in which quantum entanglement emerges 
in the study of quantum fields in an expanding spacetime. Using the tools of quantum information science, we analyse 
the appearance of quantum correlations in quantum fields due to gravitational expansion of the underlying spacetime, 
paying attention to the role played by such correlations in the current cosmological models. We discuss whether these 
quantum correlations may be present today as a remanent of the primitive universe. We also discuss the possibility 
of extracting this entanglement through local interaction with particle detectors and consider how such experiments 
could provide information about the cosmological parameters of the universe. We conclude with a discussion of possible 
applications and future research avenues in this direction. Natural units are used throughout: h = G = c = fcg = 1. 

2 Introduction to particle creation in expanding universes 

Let us consider a universe undergoing a homogeneous and isotropic expansion. This kind of spacetime is very well 
described in terms of the well-known Friedmann-Lcmaitrc-Robertson- Walker (FLRW) metric, 



where dT, is an element of a 3-dimensional space of uniform curvature, either elliptical, hyperbolic, or Euclidean. 
Expanding dY, gives 



where dQ is the line element in the unit sphere and k characterises the curvature of the space (negative, zero, or positive 
for hyperbolic, Euclidean, or elliptical, respectively). In the case where the spatial geometry is Euclidean (k = 0) we 
can write this metric as 



Even if the spatial geometry of the universe is not exactly Euclidean, the spatial curvature is close enough to zero to place 
the radius at approximately the horizon of the observable universe or beyond, so this consideration seems reasonable, 
and we will admit it from now on. 



ds 2 = dt 2 - [a{t)] 2 dZ 2 , 



(1) 




(2) 



ds 2 = dt 2 - [a(t)] 2 (dx 2 + dy 2 + dz 2 ) . 



(3) 
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2.1 Scalar field quantisation in general spacetimes 

Before considering quantum fields living in arbitrary curved spacetimes, let us quickly review for the reader the well- 
known process of field quantisation in flat spacetime. This will serve a double purpose: to introduce the notation that we 
will use in following sections and to highlight some aspects of the process of field quantisation that become relevant when 
we apply it to arbitrary backgrounds. Let us consider an incrtial observer in flat spacetime whose proper coordinates are 
the global Minkowski coordinates (t,x,y,z). He wants to build a quantum field theory for a free masless scalar field $ 
whose equation of motion is 

(□ + ?77 2 )$ = . (4) 

Since Minkowski spacetime admits a global timelike Killing vector dt, the solutions to this equation can be classified 
into positive and negative frequency. A positive frequency solution of Eq. (j4]), Uk{x,t), therefore satisfies 

d t u k (x,t) = -iujkUk(x,t) . (5) 

Importantly, this criterion would be the same if, instead of t, we consider the proper time of any inertial observer. 
Treating $ as a classical field for the moment, we can expand an arbitrary solution $>(x,t) of Eq. ((4]) as a sum of 
these positive-frequency and negative-frequency solutions with c-number coefficients. Hence, we can express <£>(a:,r.) as 
a combination of positive-frequency m(x, t) and negative-frequency u*(x,t) solutions of Eq. ((4|): 

$(x,t) = ^2[a t u t {x,t) +a*u*(x,t)} . (6) 

i 

This definition will agree with the definition for any arbitrary inertial observer Q 

The solutions Ui{x,t) (and their complex conjugates) form an orthonormal basis of solutions with respect to the 
Klein-Gordon scalar product, defined through the continuity equation to be 

(uj,Uk) = -i / d 3 x (ujdtul - u%d t Uj) , (7) 

which is positive definite over the space of positive energy solutions and satisfies 

(ui,u 2 ) = (u 2 ,ui)* , (u*,ul) = -(«2,ui) • (8) 
The modes Uj fulfil the orthonormality relations 

(u j ,u k ) = 5 jk = -(u* j ,u* k ), (u j ,u* k ) = 0. (9) 

We can now construct a Fock space for the corresponding quantum field by following the standard canonical field 
quantisation scheme. We promote the classical Klein-Gordon field to a quantum field operator satisfying the standard 
equal-time commutation relations 

[$(x, t),U(x', t)] = iS(x - x') , 

[*(»,*), *(*',*)] = [U(x,t),U(x',t)]=0, (10) 

where H(x,t) = dt$(x,t) is the canonical momentum associated with the field operator This means that we 
replace the complex amplitudes Oj and a* by annihilation and creation operators a, and a\ , which inherit the standard 
commutation relations [aj,aj] = (ui,Uj) = Sij, [ai,aj] — [af,at] = 0. We now construct the Fock space: first we 
characterise the vacuum state of the field as the state which is annihilated by all the a^. 

aj|0) = Vi. (11) 

Then we define the one-particle Hilbert space by applying the creation operators a\ to the vacuum state, 

[1^> = «1 [0> , (12) 
and by repeatedly applying the particle creator we build the complete Fock space: 

k 1 , , . . . , n* ) = 1 -J a] )"' (at )" 2 . . . (at )"* 10) . (13) 

yn\n\ . . . n ! 

Due to the symmetries of Minkowski spacetime, this quantisation procedure is independent of the particular choice 
of the inertial observer. Any other choice of time t is related to this one via Poincare transformations which do not modify 
what we would define as positive and negative frequency modes. As a consequence, the field expansion Eq. (jSJ) can be 
performed cquivalently for any inertial reference frame, and the splitting between positive and negative frequency modes 

f For notational convenience, we are using the sum symbol even when integration over frequencies is meant. In this case, 1— > f d D ~ 1 k, 
where d D ~ 1 k is the relevant integration measure. 
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is invariant. Hence, the vacuum state is also Poincare invariant, and the construction of the Fock space is equivalent 
for any inertial observer. For general spacetimes, however, we cannot assume that we have global Poincare symmetry, 
and thus we lose the objective way to build quantum fields. Instead, the vacuum will not be unique, and the notion of 
"particle" necessarily becomes observer dependent. 

If we continue analysing the scalar field but in a general background, we have to substitute the D'Alambcrtian in 
Eq. Q with the covariant D'Alambert operator, which we will denote as D g and which has the form 

n g = ^=d^y/\f\g^du, (14) 

so that Eq. Q now reads, quite simply, 

(□ 9 + m 2 )$ = 0. (15) 
Here, we have assumed that there is no coupling of the field with the spacetimc scalar curvature. This is called minimal 
coupling. But Eq. (|15|) is not the most general equation for a scalar field in a curved spacetime background. In fact, there 
are some interesting aspects of the scalar field equation in curved backgrounds with more general couplings (especially, 
conformal coupling) that will be dealt with in Sec. 12.41 but we will keep the equation in its simplest form for now. To 
extend the Klein-Gordon product, Eq. (|7|), to curved spacetime, we need a complete set of initial data — in other words, 
a Cauchy hypersurface E over which we can extend the integral^ 



(uj,u k ) = —i J dEn^ (ujd^ul - u* k d^Uj) , (16) 

where dE is the volume element, and is a future-directed timelike unit vector orthogonal to E. 

The operative question is whether the spacetime is stationary or not. If it is, then it has a timelike Killing vector 
field and we have a natural way to define positive frequency modes (xij) in an analogous way as we did for flat 
spacetime in Eq. ([5]): 

£ M V 'nUj = —iujjUj , (17) 
where ujj > and is the covariant derivative. Of course, we can construct a local set of coordinates whose timelike 
coordinate is the Killing time r associated to the isometry £ M such that it satisfies ^V^t = 1- For spacetime, which 
is a particular case of stationary spacetime, the role of r is played by the coordinate t. In contrast, non-stationary 
spacetimes do not admit a global timelike Killing vector field, and therefore there is no natural way to uniquely 
distinguish positive- and negative-frequency solutions of the field equation. In the absence of this metric symmetry 
there is an ambiguity when it comes to defining particle states. Without a natural way to split modes into positive and 
negative frequencies, there is no objective way to construct a unique Fock space, starting from the fact that there is 
no unique notion of a vacuum state. This means that there is no unique, observer-independent notion of the particle 
content of a quantum field on a non-stationary background. Nevertheless, wc can still find an approximate particle 
interpretation when the spacetime posses asymptotically stationary regions, as we will see in the following section. 

2.2 Bogoliubov transformations in asymptotically stationary spacetimes 

There are some scenarios in which the spacetime is not stationary but possesses stationary asymptotic regions. This is 
the case, for example, of some models of expansion of the universe |23j or the stellar collapse and formation of black 
holes [H]. 

Consider a spacetime that has asymptotically stationary regions in the past and in the future. We will denote 
them 'in' and 'out' respectively. In these regions we can give a particle interpretation to the solutions of the field 
equations. Namely, we can build two different sets of solutions of the field, the first one {u 1 ? } made of modes with 

kj 

positive frequency uij with respect to the comoving time in the asymptotic past. The second set {u k nt } consisting of 
modes with positive frequency u>j with respect to the comoving time in the future. In this fashion we can now expand 
the quantum field in terms of both sets of solutions to the field equation 

" (%,in< + a l it in U t) = E + 4<,out«Sr) • (!«) 



$ - 



What is more, since both set of modes are complete, one can also expand one set of modes in terms of the other by 
means of the Klein-Gordon scalar product: 

E , in out\ out / in out*\ out* /1 n\ 

[{u kt ,u k] )u k] -{u- ki ,u k . )u k . j , (19) 



C = E 



[(u k f,ul^ul-(u k fypuf*]. (20) 

Note that Gauss's theorem guarantees that the product is independent of the choice of the Cauchy hypersurface E. 
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We define the Bogoliubov coefficients as 

Then, using the properties of the Klein-Gordon product, Eq. (J5J), we can rewrite Eq. (|19p as 

„,in \ ^ I* „,out a „,out* /r,o\ 

M fe, = 2^ [ a ii M fc, ~ J ' (22) 

3 

^ = Eh'<+M^] ■ (2 3 ) 

Now, we can expand the particle operators associated with one basis in terms of operators of the other basis 

a fc„in = E { n ■'">■■ + ^< fc,,out) ' (2 4 ) 

3 

a ki ,ont = E ("-.."A: - ;,,) • (25) 

3 

Now let us consider the vacuum state in the asymptotic past region |0) in , which is annihilated by ■ V&j. We 
would like to know the form of the state |0) in in the basis of solutions of the Klein-Gordon equation in the asymptotic 
future. To compute this we use the fact that i |0) in = 0. If wc substitute . in terms of 'out' operators using 
Eq. (|2"3|). we obtain 

E(^^,out+/3;,4 jjOUt )|0) in = 0. (26) 

3 

Without loss of generality, we now make an ansatz for the general form of the state |0) in in terms of the 'out' Fock basis 
as a sum of its n-particle amplitudes 

|0) in = C|0) out + \*) ji+ C™> |*> ilA + ... + C' w » l*>, Wn +••• , (27) 
where each summand has the form 

CjU -' jn \^) h ,..., jn = E ^' 1 "" J "<- 1 ,ouf-<„,out|0)ouf (28) 

Substituting Eq. (|2~T)) in Eq. (|26[) . we obtain an infinite number of constraints for the coefficients C^i' - '- 7 ™ . We can 
readily see that the zero-particle component of Eq. (|26| can only be obtained from the annihilator acting on the 1- 
particlc component of the state defined in Eq. (|27[) . thus giving the condition 

ajiC j =0 C j = . (29) 

3 

Now, the ?i-particle component (n ^ 0) of Eq. (|26|) is obtained by acting with the annihilator on the (n + l)-particle 
component of Eq. (|2"T|) and with the creator on the (n — l)-particle component of Eq. (|27p . Thus, we know that the 
coefficients C J1 ''"' Jn+1 can be written as functions of the coefficients C ,J1 ''"' J " _1 , providing a recurrence relation that can 
be used to write all the coefficients as functions of C (the vacuum coefficient). 

This, together with Eq. (|29[) . means that the 'in' vacuum evolves to a state in the asymptotic future with no odd 
number of particles components, and that the coefficients of even components in Eq. (|27[) arc related pairwise. To find 
the form of this coefficients we use the inverse of the Bogoliubov coefficient matrix ay, which is a matrix whose ij 
component is denoted (a" 1 )^-. Therefore, given arbitrary Cj and Dj, the following identity holds: 

3 ij 

After some basic but lengthy algebra, we obtain the following expression for the 'in' vacuum in terms of 'out' modes: 



|0> in = Ccxp j 4E^( a_ % a Lout<,out I 



|0) out ■ (31) 



This is a Gaussian state, with C obtained by imposing normalisation. Notice that this state is, in general, not separable. 
This means that, depending on otij and /3y, the final state of the field can possess nontrivial quantum correlations. 

If we now compute the expectation value of the number operator in the asymptotic future N£ ut = a\, out afe-, ut 
when the state of the field is the vacuum in the asymptotic past, we find that 

(N° k f)in = in(0| A^ ut |0). n = I Ail 2 • (32) 

i 

This implies that if flij is different from zero, one would observe particle production as a consequence of the expansion. 
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2.3 FLRW Universe I: Particle production in a fast expansion period 

We have discussed in Sec. 12.21 that when a non-stationary spacetime possesses asymptotically stationary regions, it is 
possible to relate field modes in the asymptotic past to modes in the asymptotic future, being able to compute the field 
quanta creation due to the gravitational interaction. This is the case for the FLRW spacetime, Eq. ([3]), if we impose the 
additional conditions 

a(+oo) — > const. , a(— oo) — > const. , (33) 

with the former constant larger than the latter. Although measurements of distant supcrnovae have demonstrated that 
our universe is undergoing an accelerating expansion, to good degree of approximation, certain models fulfilling Eq. ([33]) 
(as the one that we will study below) can very well describe some inflationary scenarios. Roughly speaking, these include 
a period of slow expansion, then a period of very fast expansion, and then again a period of slow expansion. 

As can be seen from Eq. (|25|) . computing the analytic form of the particle operators in the asymptotic future in 
terms of the operators in the asymptotic past is not straightforward. There arc, however, some particular solvable toy 
models that allow for the study of fundamental phenomena. These include scalar models [25] and spin-^ fermionic 
models [23]. To illustrate this, let us study the scalar field case within an exactly solvable model in 1+1 dimensions. 

Let us consider the FLRW metric Eq. (|3J) and rewrite it in terms of the conformal time coordinate 



dr 



10 a{r) 
In doing this, we obtain 



(34) 



(35) 

[a(r/)] 2 . Let us assume the following specific form for the conformal 



ds 2 = [a( V )] 2 {dTf -dx 2 ). 

Following the standard notation, we define C(rj) 
factor: 

C(rj) = l + etanh(>77), (36) 

where e and p are positive real parameters controlling the total volume and rapidity of the expansion, respectively 
Imposing this specific form for the conformal factor makes it simple to study analytically while at the same time 
allowing us to study the fundamental behaviour as a function of the rapidity of the expansion and its total volume. 

It is interesting to note that for large e the FLRW metric with this conformal factor behaves, for < t <C p , like 
the radiation-dominated Friedman universe, with an exponentially fast approach to asymptotic flatness for t <C p ■ 
As mentioned above, we can consider this as a convenient way to approximate a spacetime that is asymptotically flat 
in the past and future but undergoes rapid inflation in between (with smooth transitions). The asymptotic flatness at 
each end allows us to meaningfully define different but physically meaningful particle states for the asymptotic past and 
future. Note that this example is analytically solvable only for the scalar field. The same C(r]) does not provide an 
analytically solvable model in the fermionic case, which is discussed in Sec. 13.21 

The Klein-Gordon equation in the asymptotic past and the asymptotic future in this particular case admits solutions 
of the form 



X m (v) e x P 
out (77)exp 



i(kx — uj rj) 



ln[2 cosh(p^)] 



A 



i(kx — ui.n) ln[2 cosh(p?7)] 

P 



where 



1 



2^ 



-.F 



X 



2-y/TrcJout 



P P 



P 2 

IWout 1 



where the F\ are hypergeometric functions, and 



Wout/in = 

It can be checked that indeed 



m 2 (l ± e) . 



2 ( w out i OJi, 



tanh (/w/)] 
— tanh(prj)] 

)• 



1 



1 



(37) 
(38) 

(39) 
(40) 

(41) 

(42) 
(43) 
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Using the properties of the hypergeometric functions, the computation of the Bogoliubov coefficients matrices take a 
diagonal form in this case ctkj = a%^3i Pkj = —fik&kj, so that in Eq. (|22|) the Bogoliubov transformations only mix 
modes of the same k: 

uf(x, r,) = a k u^(x, ry) + p k uT k {x, t?) , (44) 

where 

j uj out r([l - (i^ in /p)])r(-ia; ou t/p) 

k V c in r([i - (iu, + /p)])r(-i W+ / P ) ' 1 j 

fl = / ^out r([l - (io; in / j 3)])r(i^ ou t//j) 

Pfc V r([i + (i W _/p)])r(iw_/ P ) ■ 1 > 

It is important to be clear why the Bogoliubov transformation only mixes modes of the same momentum. This fact 
is not related to the particular expansion model choice in Eq. p6[) , nor to the dimension of the spacetime considered. 
Rather, it is a consequence only of the conformal symmetry (up to mass terms in the minimal coupling scenario) of 
the theory and the conformal flatness of the metric. This conformal equivalence relates the equation of motion in the 
FLRW scenario to its flat-spacetime form. This means that the general form of the Bogoliubov transformation Eq. (|44|) 
is valid for any FLRW universe with sufficiently well-behaved scaling factor and minimally coupled field, where Eq. (1331) 
is fulfilled [23]. 

From Eq. Q32[) we see that the mean number of particles in the asymptotic future when the field was prepared in 
the vacuum state in the past is 

(AOin = \^k\ 2 , (47) 

That can be interpreted as the creation of quanta in the mode k of the field as a consequence of the spacetime expansion. 
Notice that when m = 0, uj- m = w ut, which means that \j3 k \ = 0, and no particles are present in the asymptotic future. 
This is because, in the particular case of two spacetime dimensions, for m = 0, the theory is also conformally invariant. 
Let us study this case in more detail. 

2.4 FLRW Universe II: Particle detection in the conformal vacuum 

So far we have only considered the case where the field is not coupled to the scalar curvature (Ricci scalar) of the 
spacetime. This is known as minimal coupling to curvature. Indeed, Eq. (|I 5|) comes from a particular case of the more 
general action associated with a scalar field in a curved spacetime background of dimension D, which is 



S = j d D x^- [g^dpSdu* - (m 2 + £i?)$ 2 ] , (48) 

whose associated equation of motion is 

(□ B + m 2 + £R) $ = , (49) 

where £ is the coupling of the field to the Ricci curvature scalar R. In the section above, we have considered the case 
where £ = (a.k.a. minimal coupling). However, there is another case of great interest in cosmological scenarios given 
that the FLRW metrics are conformally flat. This case is that of conformal coupling. 
Consider what happens if we make a conformal transformation of the metric: 

g'^ = Mx)] 2 g^ . (50) 
The geometrical quantities in Eq. (|48|) are then changed as follows: 

g"* v =[Q(x)]- 2 g^, (51) 

VW\ = Mx)] D y/\i\, (52) 

R' =n- 2 [R-2(D -l)a g \nn + (1 - D)(D -2)g^{\nil)^{\nil)^} . (53) 

Through simple dimensional analysis we find that the field has dimensions of [length] and hence it is also rescaled 
under a conformal transformation: 

= fi(2-D)/2 $- (54) 

Now, if we want the action in Eq. (|48p to be invariant under this conformal transformation (except for a boundary 
term), we find that there are three contributions spoiling such invariance: the mass term, the derivatives in the kinetic 
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term (if fi is not constant), and the coupling to the curvature term. If we look at the massless case, however, we can 
always choose the coupling to the curvature such that the action becomes conformally invariant — namely, 
_ (D-2) 

This case (with m = 0) is called conformal coupling. In two spacetime dimensions, we see that if to = the minimal 
coupling £ = gives conformal invariance of the action, so it coincides exactly with the conformal coupling. In contrast, 
in four spacetime dimensions, we require to = and £ = 1/6 to have a conformaly invariant action. This is relevant 
to our interests since, if we apply the conformal transformation to the field equation Eq. (|49|) with D = 4, m = 0, and 
£ = 1/6 (conformal coupling) , then we obtain 

□ s + ^)* = tt 3 (n 9 < + ^')$'. (56) 

Hence, if $' is a solution of the conformally transformed equation, then $ = 0$' will be a solution of the original 
equation. 

Now, if we consider the FLRW metric from Eq. ([3J) and rewrite it in terms of the conformal time coordinate 

V=f J in> (57) 
Jo a i T ) 

we obtain 

ds 2 = [a{ri)] 2 {dn 2 - dx 2 - dy 2 - dz 2 ) . (58) 

A conformal transformation f2 = [a(?y)] _1 transforms this metric into the Minkowski metric so that g' = ij^. For rj^ v , 
R' = and the field equation is just the usual wave equation, which admits solutions of the well-known form 



Vk oc exp[i(fc • x — = exp 



k ■ x — J oj(t) dr 



(59) 



where we have defined u(t) — \k\/a(t). Using Eq. (|56|) . we have that the original equation admits solutions of the 
form Uk = [a{rf)]~~ 1 Vk- Now, Uk are positive frequency solutions at early times, so the field admits the expansion 



a: . ut, +al u* t ) . (60) 



But Uk are also positive-frequency solutions at late times, so the field also admits the expansion 

* = [ a ki,ontU ki + 4,,out u fc s ) ■ ( 61 ) 

i 

In this case, the Bogoliubov transformation is trivial, and afe ii0 ut = flfci.in, so that the early- and late-time vacua are 
the same. Indeed, for this field, the particle concept is well-defined at all times, and there is therefore a natural choice 
for the vacuum state [13 |2S]. This state is called the conformal vacuum As an obvious consequence, 

the Bogoliubov coefficients (3k = for all fc, so no particles are created in the massless conformally coupled field case. 
Therefore, modes that have a positive frequency with respect to the conformal vacuum at one given time remain so for 
all time. Hence, if the field has no particles to begin with, expansion will not create any in the future either. 

Docs this mean that, if we consider a field in the conformal vacuum, no particles at all will be detected by an 
observer? To answer this, we must consider the fact that, while the choice of the conformal vacuum for a conformally 
invariant is a natural one, it is still not unique. Even in Minkowski spacetime, the Minkowski vacuum — while unique for 
all inertial observers — is not the vacuum for an accelerating observer [5^1 [3D] . Such an observer would detect particles 
even if the field were in the "natural" choice for the vacuum state, i.e., the Minkowski vacuum |31j . This same observer- 
dependence of the vacuum persists even when the usual ambiguities associated with curvature are tamed by conformal 
invariance. As with inertial observers in fiat spacetime, however, there exists a family of privileged observers in a FLRW 
universe: those observers that see the isotropic expansion from their proper reference frame. These are called comoving 
observers. The proper time of comoving observers does not coincide with the conformal time. This means that such 
detectors actually detect particles even in the conformal vacuum. 

A particular well-known example of this is a single particle detector in a de Sitter universe coupled to a quantum 
scalar field in the conformal vacuum [T^]. We will use this example to introduce a particle detector model. Consider 
the well-known Unruh-Dcwitt detector model [33] that describes the interaction of a two- level quantum system with 
a scalar field. In its simplest form, this is a pointlikc particle detector whose interaction Hamiltonian with a scalar 
field $(x) is 

Hj = 0(t) to(t) $[a;(r )] , (62) 
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where r is the proper time of the detector, 0(t) is some appropriate switching function, m(r) is the monopole momentum 
of the detector, and $[x(r)] is the field operator evaluated along the trajectory of the detector. The monopole momentum 
operator can be written in terms of the ladder operators <r + and a~ . We can choose to expand the field operators in 
terms of positive- and negative- frequency solutions of the Klein-Gordon equation in conformal coordinates Uk(x, rj) and 
U U X ,V), yielding 

ff/ = .ge(r)y"dfc(a+c i ^+ < 7-e- i ^)(44[x(T),ry(r)]+ afcWfc [x(r),7y(T)]) . (63) 

Now, considering a massless field conformally coupled to the curvature and prepared in the conformal vacuum, one 
can readily obtain that, at the first order in perturbation theory and for long interaction times 0(t) « 1, the response 
function of a comoving Unruh-Dewitt detector is [25] 



1 f f -iB^xAWOd,/' 

FW = ^]*,jH j^— ? , (64) 

which is nonzero in general. Hence, a comoving detector will click even if the state of the field is the conformal vacuum. 
A paradigmatic example of this is the de Sitter universe, in which a comoving detector gives a thermal response [13] . 
The characteristic temperature of the detected particle statistics is known as the Gibbons-Hawking temperature and is 
proportional to the expansion rate. We will have more to say about the Gibbons-Hawking effect in Sec. 0J 

3 Entanglement generation in bosonic and fermionic free fields during inflation 

As analysed in section [2~3l when the conformal symmetry is broken there is a net particle production due to the expansion 
of the universe. As studied in Ref. [14] for bosonic fields and further analysed in Ref. [15] for both bosonic and fermionic 
fields, the generation of particles due to the expansion of the spacetime also generates entanglement in the final states of 
the field. These quantum correlations were shown to contain information about the expansion, enabling the possibility 
of deducing cosmological parameters of the underlying spacetime from this entanglement. Although arguably difficult 
to detect in a practical cosmology experiment, this novel way of obtaining information about cosmological parameters 
could provide new insight into the early universe, both theoretically and experimentally Theoretical cosmology must 
embrace the possibility of entanglement emerging as a purely quantum effect produced by gravitational interactions 
in an expanding universe and playing a fundamental role in the thermodynamic properties of FLRW spacetimes [34j . 
Observational cosmology can now look for witnesses of purely quantum effects in the early universe. In addition to this, 
a parallel research effort has arisen that uses laboratory analogues of expanding spacetime to study entanglement. This 
is discussed further in Sec. 16.21 

Remarkably, Ref. [15] reveals strong qualitative differences between the bosonic and fermionic entanglement 
generated by such expansion. The particular way in which fermionic fields get entangled encodes more information 
about the underlying spacetime than that of the corresponding bosonic case, thereby allowing us to reconstruct the 
parameters of the history of the expansion more easily. This highlights the importance of bosonic/fcrmionic statistics to 
account for relativistic effects on the entanglement of quantum fields. This difference between fermions and bosons was 
already proven to be important in other relevant phenomena in relativistic quantum information, such as entanglement 
in non-incrtial frames (35j [36} [37 , 38, 39 , stationary black- holes [40], and stellar collapse scenarios [T6] . 

3.1 Entanglement in a scalar field in a 1+1 FLRW universe 

Considering the 1+1 conformal metric in Eq. ([35]) with conformal factor given by Eq. ([36]) . we can use the solution for 
the Bogoliubov coefficients given by Eq. ([44]) directly in the formula for the 'in' vacuum state, Eq. ([3~T]) . Doing so, we 
find that the vacuum state of the field in the asymptotic past has the form, in the asymptotic future, of a collection of 
independent two- mode squeezed states [41] of right-moving and left- moving modes with momentum k: 



|0> in = (g 

k 
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L^\ n * ) ^out \ n -k)out ■ ( 65 ) 



This means that the vacuum state in the asymptotic past evolves into a mode-wise entangled state in the asymptotic 
future. To quantify this entanglement, we can study the entanglement entropy of the state in Eq. (|65]l . Entanglement 
entropy is only defined for bipartite pure states and is simply the von Neumann entropy of the reduced state obtained 
after tracing out one of the two halves of the bipartition |42j . In this case, the bipartition is into left- and right- moving 
modes in the asymptotic future ('out' modes). 

Notice that it is simply the choice of mode functions that determines whether the state appears entangled or not: 
'in' modes are separable; 'out' modes are entangled. But the state is the same in both cases; it is merely being expressed 
in different bases. Thus, the question of whether entanglement exists in the field docs not have a unique answer. This 



Cosmological quantum entanglement 



10 



ambiguity is ubiquitous in quantum information theory: the presence or absence of entanglement in any quantum system 
can only be defined with respect to a particular choice of tensor-product structure [43] ■ In a laboratory experiment, 
this tensor-product structure is usually determined by asking what measurements can be performed locally. In the case 
considered here, it is determined by what modes are "physically natural" in each of the asymptotic regions. Such an 
ambiguity in entanglement content should not be too surprising, since we have already seen that the notion of particle 
counts is also non- unique, being determined by what detectors are "natural" for a given observer (j] 

With respect to the 'out' modes, the reduced density matrix corresponding to a right-moving mode of momentum k 
obtained after tracing out the left-moving modes from Eq. (|65|) is 



Pk 



:i-|eB| 2 )£|^B| 2 "K)out(^l , (66) 



where Ob = cik/fik- Hence, the entanglement entropy Sb = — tr(/3fc log pk) becomes 

s » = los [ i-fei- J- (67) 

where, using Eq. (|45p and Eq. (|46p. 

|^B| 2 - Sinh2 2 (7r "- /P) . (68) 
sinh {-kuj+Ip) 

It is important to note that when the field is massless, there is no production of entanglement at all. This responds 
to the conformal symmetry of the theory in this particular case and is an artifact of taking 1+1 dimensions, for which 
minimal coupling is also conformal coupling (see Sec. I2.4[) . 

If the coupling to the curvature were non-zero or if we consider the massless case in higher dimensions where the 
conformal coupling is not equivalent to the minimal coupling, then would be non-zero, and some entanglement would 
be produced even in the massless case. Although no relevant differences would arguably be obtained in the behaviour 
of entanglement with the cosmological parameters, it may be useful to see how this entanglement depends on the 
coupling to the curvature, thus giving insight about fundamental properties of the underlying field theory. Although 
the cosmological information is codified in the generated entanglement, extracting that information — provided we could 
somehow experimentally access it — is not easy in the scalar case. This will be better understood after we analyse the 
fermionic case and we compare the amount of cosmological information codified in entanglement generated in each type 
of field. 

In Ref. [14] . the ability to recover information about the cosmological parameters, given access the generated 
entanglement, was analysed in the limit of small mass m -C 2 ( o v / e. They found that, in this limit, the total volume of 
expansion can be written as a monotonically increasing function of the entropy of entanglement 

2B 2 

e*-£|*B0S B )|, (69) 
m 

where 8b = Ob(Sb) is considered as a function of the entanglement entropy, inverting Eq. (|67[) . and Ek = \/k 2 + m 2 
is the energy of the fc-momentum mode. In this very same limit, the rapidity of the expansion can be related to the 
variation of the entanglement entropy with the energy of the studied mode: 




1+|0b(S b )| : 



fd £ in[|#B(s B )i 2 ]-r (70) 

Unfortunately, this limit of small mass is precisely when the generated entanglement is almost zero due to the conformal 
symmetry of the theory. While it is true that the entanglement codifies information about the parameters of the 
expansion, it is not clear how to extract, in the general case, the information contained in that entanglement of scalar 
fields. We will revisit these results later on in section |3~31 where we will compare these results with the fermionic case, 
showing that fermionic field entanglement is much better at codifying complete information about the expansion of the 
universe. 

|| See Eq. i'A'll and the discussion above it, as well as Chapter 3 of Ref. 1121 . for more discussion and history about the concept of particles 
in curved spacetime. 
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3.2 Entanglement in Dirac fields in a in a 1+1 FLRW universe 

The quantisation of a spin-^ field in a curved background is somewhat similar to the process described for the scalar 
field in previous sections but is complicated due to the presence of spin. Let us first notate the flat-spacetime Dirac 
matrices as {7 , 7 1 , 7 2 , 7 3 }. They obey the anti-commutation relation {7°, -y 6 } = 2r/ ab , where rj ah is the usual Minkowski 
metric. To study the Dirac field in curved spacetimes, we need to introduce the vierbein field [44]. This field consists of 
an orthonormal set of four vector fields that serve as a local reference frame of the tangent Lorentzian manifold at each 
point of spacetime such that 

<r = <4v ah ■ (7i) 

The vierbein enables us to convert local Lorentz indices to general indices. With the help of the vierbein, we can write 
the Dirac matrices 7^ in a general spacetime as a function of the local gamma matrices 

l* 1 = eZT • (72) 

These curved-spacetime gamma matrices fulfil 

{-f,-f} = 2g^. (73) 

We need to compute the covariant derivative taking into account the fact that we are working with a spin-i field. 
Since we have a spinor bundle defined over the spacetime manifold, the spin connection can be expressed in terms of 
the Levi-Civita connection T v a ^ as 

uf = eld^ h + e^ b t^ . (74) 
We can express the covariant derivative of the fermion field as 

D„ = {dp + r„) , (75) 

where 

r^ = (76) 

In the minimal coupling scenario, we can write the Einstein-Hilbert action for the Dirac field as 



S = / d D x det(e) 



(77) 



2 

This action leads to the covariant form of the Dirac equation: 

(KTD^ + m)yj = {yfDp + m)ip = . (78) 

Let us denote its particle and antiparticle solutions of momentum k in the asymptotic past, respectively, as u^(k) and 
u in(k)- I n the same fashion, the particle and antiparticle solutions in the asymptotic future will be denoted as u^ ut (k) 
and "out( fc )- 

In the particular case of a FLRW universe, the vierbein takes the very simple form e afJ/ = \J C{rf)rj a ^ , where 
C(rj) = [a^)] 2 is the conformal scale factor. Again, given the conformal symmetry (up to mass terms), we can relate 
the Dirac equation in the FLRW scenario to its flat-spacetime form. Using this we find, as in the scalar case, that 
the Bogoliubov transformation between modes in the asymptotic past and future does not mix different momentum 
solutions: 

uf n (k) = a±ut t {k) + ftu^k) . (79) 

We see that a particle mode in the asymptotic past becomes a superposition of one particle mode and one antiparticle 
mode in the asymptotic future. 

We will now assume a specific form for the conformal factor as we did in Sec. l2.3l so that the form of the Bogoliubov 
coefficients can be analytically determined. The same form that allowed us to obtain analytical solutions in the bosonic 
case will not work in the fcrmionic case. Instead we would need the vierbein field (and not the metric as in the bosonic 
case) to be proportional to 1 + etanh(/w/). This means that if we choose a form for the conformal factor given by 

C(r?) = [1 + e tannery)] 2 , (80) 

which should be compared with Eq. (|36[) . the FLRW metric acquires exactly the same interesting properties as in the 
scalar scenario. It describes a process of asymptotic flatness followed by a smooth transition to fast expansion and later 
followed by a smooth transition to asymptotic flatness again at a larger scale factor. The asymptotic flatness at both 
ends allows us to define physically meaningful 'in' and 'out' particle states. The similarity between this spin-i case and 
the scalar case allows us to directly compare the qualitative results obtained in each. 
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Following a process completely analogous to what lead us to Eq. (|4"5]l and Eq. (|4"6"]) . we obtain the form of the 
Bogoliubov coefficients in the Dirac case: 

^~ r(l - icj in /p)T(-iu} oxlt /p) 



f*t = 



Wi n r(l — iw+/p ± \me/ p)T{— iw+/p + ime/p) ' 

r([i - (iw in //j)])r(iw out / j0 ) 



^out 



(81) 
(82) 



w^ r(l + luJ -/p i ime//o)r(iw_/p + ime/p) 
In this case, the expression of the annihilation operators of particle and antiparticle in the asymptotic future is slightly 
more complicated than in the scalar case. To help, we can use the fact that the coefficients and j3t are not 
independent from and /37: 
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±. The polarisation tensor then takes the form 
X aa ,(k) = -2p 2 out K 2 out U out (-k,a')V(0,<T) 
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Min/out 

j-'oiit(fc) c) is a spinor solution of the Dirac equation in the curved background, and V(k, a) is a flat-spacetime antiparticle 
spinor [23] , 

If we operate as we did with the scalar case and consider 1+1 dimensions, the Dirac field has no freedom in its 
spin and behaves as a Grasmann field [45]. In this case, Eq. ([84]) and Eq. ([85]) take much simpler expressions, and the 
same applies to the 'in' operators in terms of the 'out' ones |15j . Then, with the same technique we used to obtain 
Eq. pip, we obtain the form of the vacuum state in the asymptotic past in terms of particle and antiparticle modes 
in the asymptotic future demanding that the operators a ln ^^ and &i n ,fc,o- annihilate the 'in' vacuum in the 'out' basis. 
This gives the following expression for the 'in' vacuum state in terms of 'out' modes: 
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(88) 



(89) 



Notice that due to the Pauli exclusion principle, each mode can only be populated by a single excitation (mathematically 
implemented in the anticommutation relations between fermionic operators). Following identical procedures as those 
used in Sec. 13. 1[ the entropy of entanglement was obtained in Ref. [15] as 



S F = log 
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where 1 6> F | 2 can be simplified, using Eq. flST} and Eq. ([52]), to 
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(90) 



(91) 



As in the scalar case, we again see that the entropy is zero when the mass of the field vanishes, this being the result 
of an equivalence between minimal and conformal coupling in 1+1 dimensions for Dirac fields, as well as bosonic scalar 
fields. 
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3.3 Statistics and Information encoded in entanglement 

We have seen that both scalar and Dirac field modes become entangled due to the expansion of the universe when the 
conformal symmetry is broken. In both cases, the amount of entanglement generated in the field due to the expansion 
codifies information about the underlying spacetime. However, as analysed in Ref. |15j . there are fundamental differences 
in the spectral behaviour of that entanglement. In Fig. [TJ it is shown how the entanglement depends on the momentum 
of the field mode for different values of the parameter p. For the Dirac case, entanglement peaks at a certain momentum, 
while for the scalar field, entanglement monotonically decreases with momentum and has its maximum at |fc| = 0. In 
other words, in contrast to the scalar field, for the Dirac field, there is a privileged value of \k\ for which the expansion 
of the spacetime generates a large amount of entanglement. Modes of that characteristic frequency are far more prone 
to entanglement than any others. 

In Ref. |15j . an attempt to interpret this phenomenon is given on the basis that the optimal value of |fc| can be 
associated with a characteristic wavelength (proportional to that is increasingly correlated with a characteristic 

length of the universe. Intuitively, fermion modes with higher characteristic lengths are less sensitive to the underlying 
spacetime because the exclusion principle impedes the excitation of very long- wavelength modes — i.e., those whose 
|fc| — > 0. For bosons, where this constraint does not exist, the entanglement generation is higher when |fc| — > 0. This 
can be intuitively explained by the fact that modes of smaller |fc| are more easily excited as the spacetime expands since 
it is energetically much cheaper to excite smaller |fc| modes. 

Regardless of its explanation, this natural emergence of a privileged mode in the fermion cae is a phenomenon that is 
very sensitive to the expansion parameters and thus more efficiently encodes information about the underlying spacetime 
than in the scalar case. In Ref. |15) . it is observed that, for a given mass, the frequency for which entanglement peaks 
is extremely sensitive to the rapidity of the expansion of the universe, while at the same time it is almost insensitive 
to variations on the total volume of the universe. In contrast, the amount of entanglement in the maximally entangled 
mode of the field only depends on the total volume of expansion, being insensitive to the expansion rapidity. These 
two features, which are a direct consequence of the peaked behaviour of ( | fe | , m), are telling us that the information 
about all the parameters of the expansion is codified in this peak, and if we were able to probe the field and find the 
most entangled mode, this will tell us much about the characteristics of the spacetime background for a fixed expansion 
model. 

While the expansion of the universe without conformal symmetry has been proven to generate entanglement 
regardless of the nature of the quantum field, we have seen that Dirac fields seem to codify more information about the 
underlying expansion, revealing that field statistics can play a key role in the way in which the expansion of the universe 
generates entanglement in quantum fields. 

4 The conformal vacuum, particle detection, and entanglement 

We saw several examples above of the fact that the entanglement present within a quantum field is only loosely related 
to what is detected experimentally. In particular, we saw in Sec. 13.11 that in a conformally invariant setting, expansion 
produces no entanglement between the field modes and thus, supposedly, no particles. In a particular sense, as discussed 
in that section, this is correct. But in another sense — that of local particle detectors — this is incorrect since local 
detectors with a fixed resonance frequency will see a nonzero signal under expansion even under the condition of 
conformal invariance, as we saw in Sec. 12.41 

One might argue that, because there is no entanglement generated by the expansion, such detectors do not see "real" 
entanglement. But entanglement is a physical resource that can be used for quantum communication and computational 
tasks [42], so this question should have a definite empirical answer. In fact, such an answer is available: the entanglement 
is real and useful (at least in principle). In the conformal case, it is not, however, generated by the expansion; it was 
present in the field to begin with. The effect of expansion on the swapping of this preexisting field entanglement to local 
quantum systems is an interesting topic in its own right and is discussed below. 

4.1 Curvature and thermal radiation 

The fact that black holes radiate [46] due to the quantum nature of the fields that surround them came as a shock to the 
physics community, overthrowing the prevailing intuition of black holes as universal vacuum cleaners (puns intended). 
On the other hand, the result was not entirely unexpected in light of general thermodynamic arguments about black 
holes being put forth at the time [47]. Similar results followed for acceleration [30] [31] and — most relevant to this 
work — exponential expansion (i.e., a de Sitter universe) [13] . as mentioned at the end of Section [2.41 

The groundbreaking result by Gibbons and Hawking [13] showed that an inertial particle detector in an exponentially 
expanding universe responds as if it were bathed in a thermal bath (in the rest frame of the detector) with a temperature 
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Figure 1. [Based on figures in Ref. |15| ] Bosonic field (top) and fermionic field (bottom): Entropy of entanglement for 
a fixed mass m = 1 as a function of fc for different rapidities p = 1, . . . , 40, with a fixed value of e = 1. 



proportional to the expansion rate k: 

Tgh = ^- (92) 

Z7T 

Notice that this is the same formula as that for the temperature of Hawking radiation from a black hole [46], in 
which k would represent the surface gravity of the black hole. It is also the same formula as for Fulling-Davies-Unruh 
radiation [29] [30] [31] , in which k would then represent the acceleration of the detector. The common element in all three 
cases is the presence of a horizon, which hides the region behind from observation. As such, even if the overall state 
of the field is pure, it becomes mixed because of the lack of access to part of itjj| Black holes have an event horizon, 
accelerating observers have a Rindlcr horizon at a uniform proper distance behind them, and observers in an expanding 
universe perceive cosmic horizon surrounding them. 



1 This is the standard explanation, but recent results have challenged this view 1481 . 
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Returning our focus to Gibbons-Hawking radiation, it is interesting to note that the thermal form of the detected 
radiation is robust to the details of the field undergoing inflation [13] . Also of note is the fact that there is no redshift to 
the perceived thermal radiation spectrum. That is, the thermal radiation is always perceived as being in the rest frame 
of the inertial detector. This should give the reader pause when considering that two inertial detectors in relative motion 
will each perceive a thermal bath at rest in their respective frames, both being at the Gibbons-Hawking temperature. 

4.2 Entangling power of an expanding universe 

Consider the situation in which two detectors are both on comoving trajectories (which are themselves inertial) in a 
de Sitter universe. Each detector perceives a thermal state of the field in its rest frame. But such a response would also 
be obtained if the spacetime were Minkowski and the field had just been heated up to an actual thermal state in some 
rest frame, with both detectors at rest in that same frame. As such, a single inertial particle detector cannot distinguish 
an empty but exponentially expanding universe from a heated one that is not expanding at all. An interesting question 
then is whether entanglement between two such detectors can reveal a difference. 

Rcznik, Retzker, and Silman |17| showed that entanglement in the Minkowski vacuum can be swapped to local 
quantum systems by coupling such systems to the quantum field for a finite time. We call such systems "detectors" 
even though the interaction with the field is coherent. Ref. |18j builds on this idea by modifying it in two ways and 
comparing the resulting entanglement in each case. In one modification, the Minkowski vacuum of a masslcss scalar 
field is replaced with the same field in a thermal state at rest with respect to both detectors, while the metric remains 
Minkowski. In the other case, the Minkowski vacuum is replaced by a conformally invariant scalar field in the de Sitter 
conformal vacuum, with the detectors on (inertial) comoving trajectories. 

Choosing the conformal vacuum, which greatly simplifies the calculations, can be justified on physical grounds 
because it coincides with the massless limit of the adiabatic vacuum for de Sitter spacetime |12j . Thus, we can think of 
this analysis as applying to the following two ways of adiabatically modifying the Minkowski vacuum: (1) very slowly 
heating the universe to a temperature T, and (2) very slowly ramping up the de Sitter expansion rate (from zero) to 
a final value of re, with T = k/2tt, with the detectors activated only long after this smooth transition is complete. In 
both cases, each detector responds as if it were bathed in thermal radiation in its own rest frame. The temperature 
in the thermal case and the expansion rate in the de Sitter case are chosen to satisfy Eq. (|92[) . thus ensuring identical 
responses. As such, access to one detector alone cannot distinguish the two cases. 

The coupling of the detector to the field is taken to be an Unruh-Dc Witt coupling, like the one described in 
Section 12. 4| with a fixed resonance frequency f2, separation in comoving coordinates L, and time-dependent coupling 
strength taken to be a small fixed parameter times a Gaussian in (proper) time with standard deviation a <C L. While 
the detector technically remains coupled to the field for all times due to the infinite tails of the Gaussian, a smooth cutoff 
applied in the wings of the Gaussian was found not to change the results, so the full Gaussian form is used for calculational 
simplicity to approximate a coupling that is on for a time a. The goal is to find a set of parameters (£l,o~,L) such that 
the detectors (two-level systems) would become entangled through the interaction in one case (heating or expansion) 
but not the other. 

One might expect that the entanglement persists — or even grows — in the case of de Sitter expansion (since the full 
state remains pure), while it gets drowned out by actual particles in the corresponding heated case. Counterintuitively, 
the opposite occurs. The results are shown in Figure [2] There is a region of parameter space in which a heated universe 
entangles the detectors while an expanding universe at the same perceived temperature does not. As such, the two types 
of universes can be distinguished by their "entangling power" with respect to local inertial detectors, although perhaps 
not in the way one would guess. 

For it to have operational meaning, the claim that entanglement can be used to distinguish expansion from heating 
must be capable of being verified by observers within the universe. One thing to notice right away is that if the universe is 
expanding, there is a cosmic horizon distance beyond which two parties cannot communicate |13j . Distinguishing heating 
from expansion requires that the detectors be separated by more than this distance. Thus, they cannot communicate 
with each other to verify their entanglement, nor can they use the entanglement for quantum teleportation or long- 
distance quantum communication with each other, since teleportation still requires classical communication of local 
measurement results [49]. However, if we assume that the detectors are themselves mounted on satellites and that 
each satellite can communicate with a "home planet" located half-way in between them, results of local measurements 
confirming violation of a Bell inequality can always, in principle, be transmitted to a third party on the home planet 
for verification [5J3 [51] . In addition, given an initial source of entanglement between the home planet and each satellite 
separately, the entanglement extracted from the quantum field could be swapped via teleportation back to the home 
planet [3§1 [55] . This method of remote "entanglement harvesting" from a quantum field is nowhere near practical 
since the effect is unfathomably small [18], but the thought experiment shows that in principle the existence of this 
entanglement is verifiable since it can be used as a quantum information processing resource back on the home planet. 
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Figure 2. (Based on figure in Ref. I18p Entanglement profile for detector pairs in several universes — cr is detection time, 
Q is detector resonance frequency, L is detector separation. The slanted black line is the entanglement cutoff in the 
Minkowski vacuum case (entangled above, separable below). The solid red curve is the thermal Minkowski cutoff, and 
the dashed blue curve is the de Sitter vacuum cutoff, both with perceived local temperatures satisfying 2nT = 10" 3 <r~ 1 . 
The de Sitter horizon distance (10 3 <r) is given by the dotted green line. The red star indicates one particular detector 
setup that could be used to distinguish expansion from heating. 



Several results are notable here. First, in both the thermal Minkowski and de Sitter vacuum cases, the region of 
entanglement is a proper subset of that of the Minkowski vacuum case. In the thermal case, this is expected because the 
state is now more mixed, but this is somewhat surprising for the expanding case. The results in Section 13.11 however, 
gives a possible explanation: no entanglement is generated by the expansion in a conformally invariant system. While 
this is not a rigorous explanation, it accords with the intuition that the radiation detected should not be any more useful 
(in terms of generating the physical resource of entanglement) than the vacuum itself. 

The reason that entanglement is lost in the expanding case can be understood by considering the entire system 
in conformal time instead of in the detectors' proper time. Since the field is conformally invariant, the modes behave 
like ordinary Minkowski modes with respect to conformal time. Therefore, detectors with a fixed resonant frequency 
in conformal time in an expanding universe will see the same entanglement profile as that for the Minkowski vacuum 
because such detectors are tuned to a single field mode. Having a fixed frequency in conformal time necessarily means 
having a time-dependent frequency in detector proper time (for an inertial detector), and thus, a fixed frequency in 
proper time means coupling time-dependently to many different independent field modes in an expanding universe. This 
increases local noise and reduces entanglement. This method of exchanging one type of detector behavior in conformal 
time with another type in proper time was recently proposed as a way to implement the ion-trap analogue of a detector 
in an expanding universe |53j . This proposal and several other possible analogue models are compared in Section 16.21 

5 Entanglement and quantum correlations in the cosmic background 

5.1 Inflation and the seeds of cosmic structure 

The inflationary hypothesis solves numerous problems with standard cosmology, including its two original targets: the 
horizon problem and the flatness problem [T] . The horizon problem is the observation that the universe would not have 
had time to organize its own large-scale homogeneity through causal means in the time available since the Big Bang. 
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The flatness problem is the observation that the spatial curvature of the universe appears to be fine tuned such that it 
exists for ~ 10 10 years in a metastable state instead of collapsing in on itself or expanding to complete emptiness on the 
scale of the Planck time (5.4 x 10 -44 s). 

Both problems are solved by postulating that the highly disordered, volatile, "curvy and bumpy" universe 
immediately after the Big Bang undergoes a period of rapid expansion, flattening out any curvature and flinging all 
matter and energy beyond a cosmic horizon [I]. The scale factor during this period is exponential in the cosmic time 
and thus approximates a de Sitter universe. But as we know from quantum field theory, what is left after all matter 
and energy are gone — empty space — is not entirely empty. As discussed in Section de Sitter expansion will cause the 
observable portion of the universe to appear nonempty (thermal) despite the entirety of it — or at least a much larger 
portion of it — being in a vacuum state. As such, the expansion causes the universe to possess tiny fluctuations (about 
a part in 10 5 ) in an otherwise perfectly smooth background. 

These fluctuations are taken to be the seeds of cosmic structure [54]. The colloquial explanation for this is that 
quantum fluctuations in the modes of the inflaton field (the scalar field responsible for inflation) become "frozen" when 
the size of the cosmic horizon (the inverse of the Hubble parameter) becomes comparable to the wavelength of each 
mode, creating tiny fluctuations in the background gravitational field, thus providing the necessary nucleation points for 
galaxies and other cosmic structures. Therefore, inflation now also provides a means to understand structure formation. 

The main evidence for inflation comes from the cosmic microwave background (CMB), a nearly uniform thermal 
glow of the universe at a temperature of 2.725 K, first discovered by accident [55] and later studied in detail by the 
sky mapping experiment known as the Cosmic Background Explorer (COBE). In addition to revealing a nearly perfect 
blackbody spectrum [56], COBE also revealed tiny fluctuations in temperature — signatures of the tiny fluctuations that 
led to structure formation in the early universe [2]. Its resolution was too low, however, to compare with specific 
inflationary models. Mapping the details of the anisotropy of the CMB was a task undertaken by the Wilkinson 
Microwave Anisotropy Probe (WMAP), which produced images with 30 times the resolution of those from COBE [57] 
and enabled constraints to be placed on inflationary models based on comparison of their predictions of structure 
formation with the reality of fluctuations in the CMB [3]. 

5.2 Cosmic fluctuations: the quantum-to-classical transition 

While the inflationary hypothesis, broadly speaking, solves the horizon and flatness problems, its proposal for solving 
the problem of structure formation leaves us wanting for details: what specifically is the physical mechanism by which 
quantum fluctuations in a scalar field are "frozen" by inflation into eventual classical density fluctuations in the matter 
distribution of the post-inflationary universe (and, eventually, temperature variations in the CMB)? We would like to 
have a better understanding of this process using the rigorous tools of quantum information theory. 

Several avenues for pursuing this question are available |54j , and the answers given depend on the model in question 
and possibly even on one's interpretation of quantum mechanics. Some authors [58j go so far as to postulate a need 
for new physics to explain the transition, with an appeal to dynamic collapse models — in which gravity or some other 
principle induces a real, unambiguous, dynamical collapse of a superposition into one of its constituent states. We don't 
need to go this far, however, to say something meaningful about the quantum-to-classical transition needed to account 
for cosmic structure formation. By studying the entanglement between expanding modes or regions of space and by 
considering the effects of decoherence caused by interactions within the inflaton field and with other fields, we can begin 
to discuss this physical process meaningfully. 

The approach by Nambu [20] is firmly rooted in the methods of Gaussian quantum information science and employs 
a one-dimensional lattice model of a scalar field during inflation. The lattice is separated into two finite, continguous, 
nonadjacent, nonoverlapping regions, and the entanglement between them is monitored as inflation progresses. Because 
the initial state is Gaussian and the field has no nonlinear terms, it remains Gaussian throughout the inflation, and the 
criterion of positive partial transpose for Gaussian states can be used to signal decisively the presence or absence of 
entanglement [21] . The results indicate that the two regions becomes disentangled when the size of the regions exceeds 
Hubble horizon distance. This is claimed as a necessary condition for the appearance of classicality of the fluctuations by 
virtue of entangled states failing to admit local hidden variable models for their statistics [SS] ■ The model is interesting 
because it is straightforward and simple, using only Gaussian states, and is therefore amenable to direct analysis of its 
entanglement structure. But this is also a weakness in that non-Gaussian features of inflation are not captured in this 
picture [54]. Furthermore, it focuses on the emergence of classicality within the field itself, and as we have seen above, 
properties that hold for field modes may not translate directly to what a detector (of any sort) would see. 

The effect of expansion on local detectors has been discussed in Section 14.21 but the results used a conformally 
coupled scalar field, instead of a minimally coupled field, which which is usually assumed in inflationary models. The 
results of Section 13.11 indicate that conformal coupling is quite special, and they suggest that minimal coupling will 
produce a very different result. Nambu and Ohsumi [19] address this and generalise the analysis presented in Ref. [18] 



Cosmological quantum entanglement 



18 



to a minimally coupled scalar field. This analysis involves particle detectors, thus providing a stronger grounding in 
operationalism than Ref. |20j , but the basic goal of understanding is the same: to analyze the entanglement properties of 
a quantum field during inflation, cither directly in the field itself or by considering the response of local detectors. The 
results of this study for minimally coupled fields show that entanglement vanishes at the Hubble horizon distance, in 
contrast to the result for conformal coupling [TS], in which it vanishes at about twice the horizon distance. Nevertheless, 
the quantum discord J5UJ EU IE2] m the detectors — a measure of nonclassicality that may be present even without 
entanglement — remains even when separated beyond the horizon but eventually goes to zero for sufficiently separated 
detectors. These results show that conformally coupled and minimally coupled scalar fields show different behavior in 
terms of their entangling power. 

The search for signatures of inseparability in quantum fields during inflation includes a study by Campo and 
Parentani [33] . The analysis proceeds by modeling the effects of inflation on individual field modes as a two-mode 
squeezing operation |41| . Decoherence is then added to simulate the average effect of small nonlinearities in the coupling 
that would otherwise cause the evolution to deviate from a Gaussian unitarity. The results are, unsurprisingly, that the 
two-mode squeezing operation produces entanglement between modes of opposite wave vector and that under sufficiently 
low decoherence, these states will always violate some Bell inequality [55]. Furthermore, there is a level of decoherence 
above which no Bell inequality is violated. Practical detection of this violation is a near impossibility, but the study 
nevertheless shows that if decoherence induced by nonlinearities in the evolution is small enough, the quantum nature 
of the fluctuations should survive. 

The question of nonlinearities in the evolution was taken up directly by a simulation reported by Mazur and 
Hcyl [63]. Instead of just two modes, ±k, four modes are used: ±k and ±2k. Terms are added to the field Lagrangian 
that induce couplings of the form ga±2ka±k a ±k and g* a± 2]i a±] i a±k, which arc nonlinear and allow particles to be 
exchanged: two particles from mode ±k for one particle in mode ±2k and vice versa. It is then assumed that coarse- 
graining of measurements only allows the larger modes (±k) to be detected, requiring the smaller modes at twice this 
wave vector to be traced out, thus generating entropy. Once again, entropy generation — and thus randomness in a 
measurement outcome — is taken as a sign of classicality. The results show that nonlinear interactions cause decoherence 
in individual modes by allowing the swapping of quanta to modes of higher or lower energy. These results show that 
reheating (post-inflationary particle production) is not the only mechanism for decoherence, as is commonly assumed. 
Nevertheless, the decoherence predicted is quite small and is not enough on its own to explain classicality for fields that 
do not undergo reheating. 

These results point to an interesting question, which is addressed by Franco and Calzetta [B3]. The question 
comes from applying scalar and tensor perturbations to a spacetime metric. Scalar perturbations in the metric at 
the end of inflation show up as temperature variations in the CMB, while tensor perturbations are only detectable in 
the CMB polarisation [65] . Therefore, if our instruments are only measuring temperature variation (not polarisation), 
then the tensor perturbations act as an environment for the scalar ones, causing decoherence. The authors use the 
Gell-Mann-Hartle histories formulation of quantum mechanics [55J to show that this is indeed the case. 

There remains no consensus on the mechanism by which quantum fluctuations become classical density fluctuations 
during (or after) inflation. While studying entanglement and decoherence in the field allows us to at least approach 
the problem, the interpretation of such results must be done carefully because it involves, at its heart, one of the most 
controversial aspects of quantum theory: the measurement problem. 

5.3 Relation to the measurement problem in quantum theory 

Any discussion of the quantum-to-classical transition, whether in a table-top laboratory experiment or on cosmic scales, 
necessarily begins to wade into interpretations of quantum mechanics and metaphysics, and thus we must tread carefully 
to ensure we are actually discussing that which we believe we arc discussing. For instance, no one worries that a linear 
polarisation measurement made on a circularly polarised photon reveals a definite, classical answer even though the 
outcome that can only be predicted stochastically. Even if the photon is entangled with a distant one, we still get a 
definite, classical answer. "Collapse" is a part of everyday life in experimental physics. 

As such, one should not expect that the temperature map of the CMB obtained by COBE or WMAP would 
ever appear to be in a superposition of possible maps. If inflation predicts a certain power spectrum for temperature 
fluctuations of various sizes in the CMB, then we would expect that the fixed, definite, classical temperature map 
revealed by WMAP would be statistically typical if drawn from that distribution of theoretical possibilities. This is 
the nature of all stochastic prediction. Wc stress this because the usual story is that it is entanglement in the inflaton 
field — not in the electromagnetic field, which is the field being directly observed by the satellite experiments — that is 
responsible for the density fluctuations in the early universe, which then, in turn, become imprinted on the CMB in 
the form of measurable temperature fluctuations. It is not fruitful, therefore, to ponder why our instruments detect a 
definite temperature map instead of a superposition of such maps: once a measurement is made, superpositions between 
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the possible outcome states are destroyed, and a definite answer is obtained. Such is the case with our mapping of the 
CMB temperature. 

On the other hand, if the primordial entanglement from the inflaton field were somehow transferred to other fields 
that we have access to today, such as the electromagnetic field, then it is not impossible (at least in principle) that we 
might observe entanglement in the individual photons of the CMB coming from different spatial directions by looking 
at correlation statistics over many observations, such as is done to verify entanglement in parametric down-conversion 
experiments [67l [68] . 

This is exceedingly unlikely, however, at least for the electromagnetic field, because entanglement and quantum 
coherence are extremely fragile. We do not need to look far for things that destroy them. Dccohercncc is a constant 
nuisance in experimental quantum physics, and as such, the burden of proof rests squarely on the shoulders of those who 
wish to claim that long-range coherence could be maintained for millennia after inflation ended. Even if individual field 
modes are entangled as a result of the gravitational interaction, once they interact with anything else in the universe, 
their coherence will quickly be lost unless the fields are exceptionally weakly interacting. We fully expect there are 
a number of such dccohcrcnce mechanisms through self-coupling and coupling to other fields, including gravity. It is 
important to determine which of these processes are dominant, and it is here that the community may fruitfully focus 
its efforts. 

6 Prospects: Quantum entanglement as a tool for Cosmology 

In the hopes of stimulating further research into the connection between quantum entanglement and cosmology, we 
conclude with some discussion of further research avenues in this area. We discuss the possibility of the cosmic neutrino 
background possessing entanglement from the early universe and surviving all the way to the present due to its weak 
coupling to matter. We also discuss the importance of analogue experimental models that provide a laboratory-accessible 
testbed for cosmological models. Finally, we conclude with a brief summary of possible future prospects for using 
quantum entanglement as a tool for theoretical and experimental cosmology. 

6.1 Quantum correlations in the cosmic neutrino background 

We have seen in previous sections that the quantum correlations present in quantum fields in the early universe would be 
hard-pressed to survive to the present due to decoherence induced by field interactions. This means that if remnants of 
these quantum correlations are to survive and be detectable today, they must do so in those quantum fields that interact 
very weakly. In this regard, the CMB is highly non-optimal since the electromagnetic field interacts with everything 
that has electric charge, thus leaving a background that is almost certainty decohered by now. Alternatively, the cosmic 
neutrino background |69j does not have this problem since neutrinos rarely interact with anything at all. As such, the 
neutrino background might prove to be a better candidate for residual quantum correlations than the CMB. 

Due to the fact that these low energy neutrinos only interact weakly with the rest of the quantum fields they are 
extremely difficult to detect so that, unlike the CMB, the neutrino background is not easily directly observable. There is, 
however, indirect experimental evidence of its existence [70], and it is estimated that today the neutrino background has 
a temperature of roughly 2 K. On the other hand, the neutrinos' introverted nature, which thwarts our efforts to detect 
them, nevertheless becomes a feature when it comes to the study of quantum correlations in the primitive universe, 
since the strong quantum correlations present in the early universe would have been protected from decoherence by a 
lack of interaction with other fields. 

Despite the difficulty of its detection, if our sole criteria is to improve on the CMB as a possible source of primordial 
cosmic entanglement, then the neutrino background remains a candidate. Searching for entanglement in the CMB is 
severely hindered by the fact that the light we observe today dates back only to recombination — the formation of the 
first stable atoms and the resulting decoupling of matter and radiation that followed — which occurred approximately 
~ 10 5 years after inflation was to have ended. That's an impossibly long time for entanglement to survive under strong 
coupling to matter. Neutrinos decoupled from matter much earlier than the electromagnetic field — approximately one 
second after inflation |71] . As such, the required coherence time for any quantum entanglement to survive in the neutrino 
background (assuming no decoherence during its travel to our detectors) is 12 orders of magnitude lass than what would 
be required for entanglement to survive in the CMB (one second versus 10 5 years). Despite the fact that the idea 
of probing the neutrino background for primordial cosmic entanglement is nowhere near a practical proposal, these 
arguments suggest we should not dismiss the possibility out of hand. 

6.2 Analogue models 

Most predictions of curved-spacetime QFT, including nearly all of the effects discussed in this paper, are extremely 
difficult to test directly, confined as we are to a low-curvature region of spacetime, with dimensional analysis conspiring 
to make the most straightforward effects extremely tiny. This situation has given rise to a number of proposals for 
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testing analogue quantum field theory in curved spacetime, the seminal example being Unruh's proposal for detection 
of the sonic analogue of Hawking radiation in supersonic fluid flow |72j . 

In such experiments, spacetime is replaced by a system of coupled atoms, either as discrete entities or in a fluid, and 
field modes are replaced by the collective normal modes of oscillation of the atoms — effectively a "phonon field" [75] . 
The advantage of this over probing fundamental quantum fields in outer space is that, while the observable effects in the 
laboratory setup are conceptually the same as (analogous to) those predicted for quantum fields on a curved background, 
the parameters of the experiment may be adjusted to give much stronger and more easily observed effects. The most 
common analogue curved-spacetime field theory proposals involve Bose-Einstcin condensates (BECs) [74] [75] [76] [77] 
or ion traps [73 [53] as the analogue spacetime. We say a few words about these here and refer the reader to the 
aforementioned references, plus the review articles of Refs. (79] [73], for more details. We also note that liquid helium 
has also been proposed as an analogue testbed [50], as well as schemes based on photons in a microwave guide and 
superconducting circuits |81) . 

BECs provide a natural system for studying the analogue of FLRW expansion [74] [75] [76] [77] . The condensate 
field is linearized into a classical background plus a quantum field of small fluctuations that lives on top of it. When 
a fully condensed system is allowed to expand in the laboratory, the emulated effect is the expansion of the spacetime 
on which the quantum fluctuations live. Furthermore, as expansion proceeds in the laboratory, the initial correlations 
in the ground state of the field naturally evolve into classical fluctuations in the density of the background condensate, 
which can be observed simply by taking a picture of the condensate during expansion. In a freely expanding BEC, a 
sonic horizon is formed, allowing an experimental study of the "freezing" of quantum fluctuations into a classical density 
distribution as the size of a given mode stretches beyond the sonic horizon. 

The main advantage of trapped-ion analogue gravity experiments [78[ I53j over other proposals, such as BECs, lies 
in the exquisite quantum control over state preparation, interaction, and readout [55]. Each ion is equipped with its 
own motional detector in the form of laser-induced coupling of its vibrational motion to its electronic state. Effectively, 
then, we have a discrete phonon field with a detector mounted at every point that we can turn on and off at will in a 
continuous and highly controlled manner. In addition, cooling to the ground state is achievable with high fidelity. The 
tradeoff for this precision and "digital" nature of the field and its detectors is in the small maximum number of ions 
confinable to a single linear trap (currently, about ten |82j). Two methods of using ion traps to study cosmic expansion 
have been proposed. The first |78] allows the ions to physically expand for a time, simulating expansion starting and 
then stopping, as in Sec. 12.31 Another proposal [53] uses the properties of conformally invariant expansion to encode all 
details of the scale factor in the parameters used to couple the electronic and vibrational degrees of freedom. This allows 
detector responses during expansion of a conformally coupled field to be emulated in the ion trap without physically 
moving this ions at all. (Similar techniques are used in BECs, as well [7JJ US].) In addition to these applications, a 
linear ion trap can be used to emulate entanglement swapping from a flat-spacetime quantum field to local detectors 83 . 
Ref. [53] proposed modifying this setup to perform the analogue of the thought experiment proposed in Ref. [18] and 
discussed in Sec. 14. 2\ possibly bringing entanglement-based discrimination of (analogue) gravitational parameters into 
the realm of laboratory testability. 

6.3 Other research avenues 

We end with some thoughts on further research avenues in this area. 

Coupling to curvature — An important observation is that conformal invariancc, while greatly simplifying calculations 
for expanding spacctimes, is often not present in inflationary models. Minimal coupling is used instead. While the two 
settings coincide in two dimensions, they differ in every other dimension, including that of our own 4-dimcnsional 
universe. There are substantial differences between the physics of the two types of coupling, even in the massless case, 
as we have stressed many times throughout this article, with the most important being that expansion of a conformally 
coupled field generates no entanglement in the field modes and no particle production if the expansion ceases (Sec. 13.1]) . 
Therefore, analyses that include the ability to select the conformal coupling constant (or at least to allow it to be set 
to zero) are a natural generalisation of many of the results presented here and would allow us to examine how the 
entanglement produced scales with this coupling strength. 

Entanglement extraction — In addition to the entanglement produced in the field modes themselves, we can also 
ask how the field coupling to curvature affects entanglement extraction from the field using local detector couplings. 
Ref. [19j pursued this as a first step, generalizing the results of Ref. [18j from conformal to minimal coupling. But 
we know from analogue proposals [74[ I75[ 153] , as well as from more recent results showing that entanglement exists 
between [M] and can be extracted from [55] the past and future wedges of a spacetime (instead of from the left and 
right wedges), that modulating a detector's resonance frequency and interaction strength can be used to optimize the 
extracted entanglement. Even the early proposal for this procedure, Ref. [IT] , introduced a rather complicated time- 
varying coupling to amplify the amount of entanglement extracted from the vacuum state. In addition to optimising 
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the extracted entanglement, we might also ask what time dependence of the interaction would optimise the extraction 
of information about cosmological parameters and spacetime structure using this method. 

Anisotropy of inflation — Previous literature analysing the emergence of entanglement in quantum fields during 
inflationary periods presented results based on the assumption of isotropy in the expansion. The question remains open 
as to whether the entanglement generated during inflationary periods may be sensitive to possible anisotropics in this 
expansion. It is known that particle creation feels the difference [53], but it still remains to be explored to what extent 
this happens and how much information about the anisotropics would be encoded into the entanglement generated in 
the inflationary period. 

Analogue models — Perhaps one of the most promising avenues of research is in quantum emulation of expansion 
through analogue experiments. In the ion trap architecture, the approach of Ref. [53] is complementary to that of 
Ref. [78| in that the former modifies detector parameters to emulate expansion of a conformally coupled field, while the 
latter literally expands the trap, allowing the ions to manifest actual phonon production with expansion over a finite 
time. Perhaps the two could be combined to study more complicated phenomena. Alternatively, Ref. |87j proposes a 
BEC model in which atoms in two different hypcrfine states are trapped together and allowed to emulate more complex 
inflationary scenarios. The development of a reliable, robust, and highly tunable laboratory testbed for analogue inflation 
would be of great experimental value in testing predictions in theoretical cosmology. 

Observational astronomy — Our best data sets today regarding the early universe come from observations of the 
CMB |57j . The possibility of using the cosmic neutrino background to search for primordial entanglement was already 
discussed in Sec. 16.11 Since this proposal remains far beyond current technology, the CMB is our best hope for testing 
theories of cosmology — through the temperature maps produced by WMAP [3] and through ongoing efforts to map the 
polarisation of the CMB, which encodes cosmic information not found in the anisotropics alone [65j . 

The research effort into the connection between cosmology and entanglement is vast, ranging from studies of 
entanglement in field modes, to the ability to swap this entanglement to local detectors, to the extraction of cosmological 
parameters from such measurements, to direct comparison of observational data and theoretical models, to the 
experimental probing of theoretical proposals in analogue experiments. With multiple fronts opening into unexplored 
territory, we hope the reader will find the results illuminating and inspiring of further research into the connection 
between cosmology and entanglement at the theoretical, observational, and laboratory-experimental levels. 
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